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ABSTRACT 


Strain energy release rates were used to correlate cyclic debonding 
between metal sheets and composite reinforcement. An expression for the 
strain energy release rate was derived and applied to fatigue test results 
for three material systems : graphite bonded to aluminum with both a room 

temperature and an elevated temperature curing adhesive, and S-glass bonded 
to aluminum with an elevated temperature curing adhesive. For each material 
system, specimens of several thicknesses were tested with a range of 
fatigue loads. Cyclic debonding was monitored using a photoelastic tech- 
nique. A close correlation was found between the observed debond rates 
and the calculated strain energy release rates for each material system. 



INTRODUCTION 


Adhesive bonding is becoming widely used in aircraft structures for 
joining structural components and for making efficient materials. In joints, 
it eliminates severe stress concentrations that are usually introduced with 
mechanical fasteners. Although it can be advantageous in met al-t o -metal 
joints, it is particularly applicable to Joining composite materials, whose 
static strength is very sensitive to stress concentrations. Also, bonding 
is used to make more efficient structures by joining separate materials to 
form one system. For example, hybrid systems, formed by bonding metal and 
composite layers, have higher static strength — for equal weight and stiff- 
nesses [l] — than metals while being more reliable than composites alone. 
Indeed, composites themselves derive their high efficiency from bonded col- 
lections of constituent materials . 

When these bonded structures are subject to fatigue loading they may 
be susceptible to a little-considered mode of fatigue failure; cyclic de- 
bonding [2]. Because most practical structures will be subject to fatigue 
loading, cyclic debonding should be considered in their design. Unfor- 
tunately, a designer has virtually no rationale to account for cyclic de- 
bonding; consequently, either the reliability or efficiency of his structure 
may suffer. As a first step in supplying a design rationale, the objective 
of this paper is to present a cyclic debond analysis. 

The analysis was developed for simple laminated specimens made of 
aluminum alloy sheet bonded to graphite or fiberglass composites and tested 
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under constant-amplitude fatigue loading. It is "based on the correlation 
of the observed rate of debonding with the computed rate of strain energy 
release as the debond extended. 

SYMBOLS 

The units for the physical quantities defined in this paper are given 
in the International System of Units (Si) [3]. The measurements and calcula- 
tions were made in the U.S. Customary Units. 

a Debond length, m 

da/ dW Debond propagation rate, m/ cycle 

c,n Curve fit parameters 

2 

E Young's modulus, MN/m 

G Strain energy release rate, J/m 

L Length, m 

P Applied load, N 

R Ratio of minimum-to-maximum applied stress 

2 

S Stress in composite core of Region A, MN/m 

t Thickness, m 

At Change in temperature, K 

U Strain energy, J 

3 

V Volume , m 

w Specimen width, m 

x,y Cartesian coordinates , m 

a Thermal expansion coefficient, K ^ 
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6 Deflection, m 

e Strain 

2 

a Stress, MN/m 

3 

Strain energy density, J/m 
Subscripts : 

1 Denotes aluminum cover 

2 Denotes composite core 

A Denotes Region A 

B Denotes Region B 

C Denotes Region C 

EXPERIMENTAL PROCEDURE 
Specimens and Loading 

The specimen configuration used in the present study is shown in Figure 
1. The specimen was composed of two 7075-T6 aluminum alloy sheets bonded to 
a unidirectional composite core of graphite or S-glass (see Table i). Two 
bonding materials were used for the graphite core specimens: EPON 927 which 

cures at room temperature (material system l), and AF 126 which cures at 
39*+ K (material system 2). The AF 126 was also used to bond the S-glass core 
specimens (material system 3, see Table II). The abrupt change in section 
of the specimen purposely introduced a severe stress concentration in the 
bonding materials. Under cyclic loading, debonding started readily at this 
high stress concentration. 
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The specimens were tested axially under constant-amplitude fatigue 

loading with R = 0.1. Maximum stresses in the unreinforced composite core 
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range from 211 to 1210 MN/m . All of the tests were performed at a frequenc; 
of 10 Hz. 


Measurement of Debond Rates 

The debond front was monitored continuously by a photoelastic technique 
Photoelastic coatings were bonded to the aluminum sheets, and the specimen 
was viewed through a polarizer and quarter-wave plate. Under loading, 
isochromatic fringes developed at the debond front due to the high strain 
gradient in that vicinity. Figure 2 shows the location of the photoelastic 
coatings on the specimen and typical isochromatic fringes. The isochromatic 
fringes were photographed at specified load cycle intervals to relate the 
debond front location to the number of applied load cycles. 

STRAIN ENERGY RELEASE RATE EQUATIONS 

Fatigue crack propagation rates in metals have been correlated using 
the strain energy release rate, G, [4] where 

0 - P <£> - £ <!> 

Hie term P(d6/da) is the work done by the applied load as the crack extends, 
and dU/da is the change in strain energy as the crack extends. Intuition 
suggests that a similar correlation may be valid for debond propagation. 
Consequently, an expression for G was developed for the specimen 
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configuration used in this study. This expression was determined using a 
one-dimensional elasticity analysis described by the following discussion. 

To simplify the analysis, the specimen was separated into three regions 
A, B, and C (figure 3). In Regions A and C only uniform stresses in the 
x direction are significant. Consequently, these regions could be analyzed 
by an elementary elasticity method. Region B had a complex stress distribu- 
tion and could not be analyzed using elementary methods. The stress distri- 
bution in Region B was assumed to remain constant as the debond extended. 

As will be shown, this constancy eliminated the need to calculate the stress 
distribution in Region B. 

An expression for (ds$/da) in equation (l) was derived from the change 
in the end deflection of the specimen, , caused by an increment of de- 
bonding, da. Before debonding occurred, the end deflection was given by: 

6 = 6 A + 6 B + 6 C (2) 

and after debonding by: 

s' = s' c (3) 

The stress distribution in Region B was assumed to be the same before and 
after the region translated; thus, the 6 term did not change with de- 

-D 

bonding. To find d6, equation (2) was subtracted from equation (3) 
yielding: 

S' - S = d6 = (6^ - 6 A ) + (&' c - 6 C ) (U) 
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The deflections on the right side of equation (k) can "be expressed in the 
general form: 

6 = eL (5) 


Substitution of equation (5) into equation (U) yields: 

d6 = [e A (L A + da) - e^] + [e c (L c - da) -e^] 

d6 = (e A - e c )da 


or 


dtS 

da 



( 6 ) 

(T) 

( 8 ) 


The expression for dU/da in equation (l) was derived by calculating 
the change of strain energy in the specimen, dU, resulting from an increment 
of debonding da. Employing the same reasoning used in the development of 
equation (h), the change in strain energy is: 

U* - U = dU = (U A - U A ) + (l^ - U c ) (9) 


The strain energies on the right side of equation (9) can be expressed in 
the general form: 

U ■ <|>V (10) 
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Substitution of equation (10) into equation (9) yields: 

dU = + da) - *A2 wt 2 L A ) 

+ {[ 24 . cl wt 1 (L c - da) + <l> C2 wt 2 (L c - da)] 

- t^Ol^C + * 02 ^ 2 ^ (11) 

au = vU A2 t 2 - (* 02 t 2 ♦ 2* 0 iti)]da (12) 


or 


£ * W[ \2 t 2 - (♦ 0 2 t 2 + 


Vl 11 


(13) 


Substitutions of equations (8) and (13) into equations (l) yields: 

G “ P(E A - ^ - w[<f A2 t 2 - + 2<|, 01 t l )] (11 ° 


This equation is evaluated in the appendix in terms of applied stress, 
temperature change, material parameters, and specimen configuration and 
leads to 


t t E v 

G = V a A + W [S ' AT(ai " a2)E23 


(15) 
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RESULTS ML DISCUSSIONS 


Figure ^ shows a sample plot of debond length against the number of 
applied load cycles. The shape of this debond curve was typical for all 
three material systems tested. Initially, debonding was nonlinear with 
respect to the number of applied load cycles, but became linear as the de- 
bond progressed. The main focus of this paper is on the linear portion 
of the debond behavior; however, a brief discussion of the nonlinear region 
is merited. 

Nonlinear debonding occurred when the debond front was near the change 
in cross section of the specimen. The texture of the failure surface in- 
dicates that the failure mechanism changed as the debond extended. Figure 
5 shows a photograph of the fracture surface of a typical debond specimen 
of material system 2. The graphite-composite core is shown on the left and 
the mating aluminum cover sheet on the right. As the photograph indicates, 
the failure mechanism was initially cohesive, but changed to predominantly 
adhesive as the debond extended. Similar behavior was observed for material 
system 3. However, for material system 1 the failure mechanisms seemed to 
be reversed; initially the failure was adhesive, but changed to cohesive as 
the debond extended. For all three material systems the nonlinear portion 
of the debond behavior seemed to be related to a transition in failure 
mechanism. 

Debond propagation rates were determined from the linear portion of 
debond versus cycle plot for each specimen. Figure 6 is a plot of rate 
against G, equation (15), for the three material systems. This figure 
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shows debond rate to he a single valued function of G for all material 
thicknesses. Table 3 presents these data in tabular form. An equation 
of a form which had successfully correlated fatigue crack propagation data 
for isotropic metals [ 5 ] is: 


ff=c(o)" ( 16 ) 

This equation fits the data in figure 6 quite well (dashed curves). The 
constants c and n were determined using least-squares techniques and 
are given in the following table: 


Material System 

c 

n 




1 

1.19 x 10 y 

3.30 


-9 


2 

2,32 x 10 * 

2.15 

3 

5.63 x 10~ 9 

3.76 


Equation (l 6 ) may be useful to predict debond rates in composite reinforced 
structures subjected to fatigue loading for the systems studied. However, 
for other materials systems, appropriate values of c and n must be 
obtained from cyclic test data. 

Figure 6 shows that c and n differ for each material system. This 

difference was probably due to the different bond strength and residual 

thermal stress of each system. The residual ther mal stress in the graphite 

aluminum specimens bonded at elevated temperature was calculated to be as 
2 

high as 120 MN/m in the graphite core. Because residual thermal stresses 
of this magnitude are significant in comparison to the applied stress in 
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the core, they were included in the derivation of the strain energy release 
rate. However, "because both the residual thermal stresses and bond 
strength differed from system to system, isolation of either effect was not 
possible for the limited tests reported here. 

CONCLUDING REMARKS 

A fatigue analysis method was developed for cyclic debonding of 
laminates composed of metal and composite layers. The strain energy re- 
lease rate, G, correlated the debond propagation rates from each of three 
series of laboratory tests for several applied loads and thicknesses. Speci- 
mens tested were graphite bonded to aluminum at room temperature, graphite 
bonded to aluminum at elevated temperature, and S-glass bonded to aluminum 
at elevated temperature. The test data were well represented by an equa- 
tion of the form da/dN = c(G) n . For the systems studied, a closed-form 
expression was developed for G, which included residual thermal stresses. 

This study indicates that G promises to be a tractable tool with 
which to analyze the relations between cyclic debonding and fatigue loading 
in laminated structures . For simplicity in establishing G as a correlating 
parameter, tests were restricted to constant amplitude loading at room tem- 
perature with a simple geometric configuration. The correlation achieved 
herein was on specific material systems that may be applicable to aircraft 
structures. However, to be useful for structural applications the model 
should be refined and verified for more complex fatigue environments and 
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for complicated configurations. Also, since high residual thermal stresses 
may occur for some material systems, their effect on cyclic debonding should 
be established. 
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APPENDIX 


Development of the Strain Energy Release 
Rate Equation 

The strain energy release rate for debonding of a metal overlayed with 
a composite (figure 3) was given in the body of this paper by equation (l4): 

0 * Pte A - e C> - wt ' t ’A2 t 2 - ( *C2*2 + Vl 11 (A1) 

The strains in equation (Al) were found by requiring equilibrium, 
strain compatability, and a constitutive relation between stress and strain 
in regions A and C. Equilibrium is satisfied by: 

P = w/crdy (A2) 

For the one dimensional analysis used in this problem, strain compatability 
was assured by assuming that the strain through the specimen thickness is 
constant . The constitutive relation for the problem is given as [6] 

a = E(e - aAT) (A3) 

Using these three relationships, the strains in region A and C can be cal- 
culated. 

In region A, compatability is satisfied. Equilibrium is satisfied 


when: 



The constitutive equation (A3) for this region is; 


a A2 = V £ a - V® 1 


(A5) 


Substituting equation (A5) into equation (A^) and solving for e. yields: 


e A * * V T 


(a6) 


For region C, compatability is satisfied vhen: 


£ C1 = £ C2 £ C 


(A7) 


Equilibrium is satisfied when: 


P = St 2 v = 2t 1 wcr cl + cr C2 t 2 w 


(A8) 


The constitutive relationships for the metal and composite respectively 


are: 


°C1 * V e o - °1 4T) 


(A9) 


°C2 = V e C - “2 aT) 


(A10) 


Substituting equation (A9) and (A10) into equation (A8) and solving for 
e c lead to: 


St 2 + AT(2K L t 1 a 1 + Bgt^g) 
£ C = 2E 1 t 1 + E 2 t 2 


(All) 



The strain energy density can be expressed as: 


_ cr 2 _ TEfe - QtAT) V 


2E 


2E 


or 


(A12) 


| [e 2 - 2e«AT + ot 2 (AT) 2 ] 


For region A, substituting equation (A6) into equation (A12) yields: 

^2 = ~2 [e A " 2£ A a 2 AT + a 2 (AT)2] (A13) 

Similarly for region C, substitution of equation (All) into equation (A12) 
for the metal yields: 

E 

<t> C l = -7T [e 2 - 2E c a i AT + a 2 (AT) 2 ] (AlU) 


and for the composite yields: 




2e c a 2 AT + ct 2 (At) 2 ] 


(A15) 


Substitution of equations (A 6), (A8) , (All), (A13) , (Al^), and (A15) 
into equation (Al) yields: 


0 • t [S - 4T( “l - W* 


(AX6) 


the expression for strain energy release rate. 



TABLE I. - MATERIAL PROPERTIES 


MATERIAL 

E 

a 

MN/m 2 

K’ 1 

7075-T6 

ALUMINUM ALLOY 

71 x 10 3 

22.5 x 10" 6 

GRAPHITE -EPOXY 


EEH51 

S-GLASS-EPOXY 

61 x 10 3 



TABLE 2. - SPECIMEN MATERIAL SYSTEMS 


MATERIAL 

SYSTEM 

COMPOSITE 

BONDING 

MATERIAL 

BONDING 

MATERIAL 

THICKNESS 

CURE 

TEMPERATURE 

mm 

K 

1 

GRAPHITE -EPOXY 

EPON 927 

0.13 

RT' 1 * 

2 

GRAPHITE -EPOXY 

AF 126 

0.13 

394 

3 

S-GLASS-EPOXY 

AF 126 


394 


(1) RT = 294 K 







































TABLE 3. - STRAIN ENERGY RELEASE RATES 
AND CYCLIC DEBOND RATES 


MATERIAL 

SYSTEM 

TYPE 

THICI 

IT 

(NESS, 

m 



*2 




S (1) 

504 

706 

807 

1210 





.51 

0.81 

da/dN (2) 

0.000714 

0.092 

0.146 

2.500 







G (3) 

16.1 

31.5 

41.2 

9Z6 







S 

426 

475 

527 

590 

639 

698 

746 

TYPE 1 

1.02 

1.64 

da/dN 

0.0168 

0.0229 

0.166 

0.373 

0.498 

1.32 

1.78 




G 

23.6 

29.3 

36.2 

45.3 

53.2 

63.4 

72.4 




S 

2ii 

270 

324 

376 

429 

543 

598 


1.60 

2.50 

da/dN 

0.00689 

0.0315 

0.0569 

0.0785 

0.163 

1.32 

5.59 




G 

8.9 

14.7 

20.8 

28.0 

36.6 

58.6 

70.9 


(1) STRESS, MN/m 2 

(2) DEBOND RATE. \i m/cycle 

(3) STRAIN ENERGY RELEASE RATE 


WITH THERMAL RESIDUAL STRESSES, G, 
joules/ meter 



TABLE 3. - CONCLUDED 


MATERIAL 

SYSTEM 

TYPE 

i 

K- 

s, 

m 












*2 


— 





S (1) 

390 

445 

541 



667 1 

667 1 

mm 

716 


1.02 

1.57 

da/dN (2) 

0.015 





EilMU 




TYPE 2 



G (3) 

60 


89 



117 

117 

126 

MM 



S 


423 

478 


535 

589 

631 




1.60 

2.54 

da/dN 




0.0627 



unim 






G 

78 


121 

125 

139 

158 

173 





1 

S 

334 

390 

446 


558 


670 

725 

781 


1.02 

.86 

da/dN 

0.0305 

0.0607 

0.0772 

0.184 

0.409 

0.605 

1.26 

2.03 

295 




G 

53 

67 

83 

102 

119 

139 

162 

185 

211 

TYPE 3 



s 

217 

274 

328 

384 

438 

492 

547 

595 

658 


1.60 

1.40 

da/dN 

0.0078 

0.0254 

0.0838 

0.181 

0.546 

1.14 

1.85 

3.33 

5.92 




G 

47 

64 

84 

105 

129 

155 

184 

212 

251 


(1) STRESS, MN/m 2 

(2) DEBOND RATI, (i m/cycle 

(3) STRAIN ENERGY RELEASE RATE 


WITH THERMAL RESIDUAL STRESSES, G, 
joules/ meter 
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FIGURE 1 -SPECIMEN CONFIGURATION 
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[FIGURE 3 - TEST SPECIMEN BEFORE AND AFTER AN INCREMENT OF DEBONDING 
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FIGURE 4 -TYPICAL VARIATION OF DEBOND LENGTH WITH CYCLES 
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